Reconstruction of the environmental correlation function from single emitter photon 

statistics: a non-Markovian approach 
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We consider the two-level system approximation of a single emitter driven by a continuous laser 
pump and simultaneously coupled to the electromagnetic vacuum and to a thermal reservoir beyond 
the Markovian approximation. We discuss the connection between a rigorous microscopic theory 
and the phenomenological spectral diffusion approach, used to model the interaction of the emitter 
with the thermal bath, and obtained analytic expressions relating the thermal correlation function 
to the single emitter photon statistics. 
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I INTRODUCTION 

Single Molecule Spectroscopy (SMS) is a powerful ex- 
perimental tool with widespread applications in nano- 
tcchnology, nano-biology, quantum communication and 
quantum computation [l|-|4| . Given the possibility to iso- 
late a molecule, a quantum dot, or an atom allows for 
investigation of the quantum dynamics of the system at 
the microscopic level. The simplest method of the in- 
vestigation consists of exciting the single emitter by a 
laser and analyzing the outcome radiation. The statistics 
of the photons, spontaneously emitted by such a single 
quantum emitter, depends on its interaction with the en- 
vironment and permits extracting information on the lat- 
ter at the level of atomic distance and time scales. Thus, 
SMS is a key for nano-technology advancement [EHl] . 

The theory of an open system, extensively developed 
over the last thirty years, combines phenomenological 
and microscopic approaches, where the environment 
is typically modeled as a thermal bath of harmonic 
oscillators (ill- 21 1. An important result of the theory is 
the derivation of a reduced Liouvillc equation, obtained 
by tracing out the irrelevant degrees of freedom. The 
most general integro-differential form of this equation 
shows that all microscopic details of the coupling 
to a reservoir are unified within the memory kernel 
given by the environmental correlation function /3t(£) 
(defined below). In the so-called Markovian limit, the 
environmental correlation function may be considered 
as a (5-function of time (in comparison to the time-scales 
of the unperturbed particle, determined by the inverse 
of the eigneenergies of the system), which leads to an 
irreversible, pure semigroup evolution. In practice, 
however, the Markovian dynamics is only an approxi- 
mation, whose validity depends on several factors such 
as the density of environmental states, temperature 
and the explicit form of the coupling 19rl21[. At very 
high temperatures the inaccuracy of the Markovian 
approximation is negligible 0, [13]. However, at low 



temperatures, fast changing or vanishing density of 
environmental states, as known to occur in Photonic 
Band Gap materials, the evolution of the open system 
corresponds to the non-Markovian regime, where the 
revival effects are enhanced while dissipation effects are 
minimized 16 - 22j ■ 



Up to now, the theoretical investigations of SMS were 
mainly based on a phenomenological approach assuming 
that the coupling to a thermal bath may be imitated 
by a real noise r\t artificially added to the unperturbed 
frequencies of the emitter fi^l - EHl ] . The properties 
of this noise, called a spectral diffusion process, and 
especially its two-times correlation function (rjtT] T }, are 
meant to reflect all possible effects of the interaction 
with the surroundings. Such an intuitively pleasing 
phenomenology simplifies the analysis, however, it leaves 
unclarified the connection of the model to the rigorous 
microscopic theory and does not give access to the 
temperature dependent parameters of the reservoir. 

In this article we provide a quantitative connection 
between the spectral diffusion model and the micro- 
scopic method, thus allowing SMS techniques to provide 
detailed information on the environmental dynamics 
beyond the Markovian limit. In section II we discuss the 
theoretical limitations of the spectral diffusion approach 
through the comparison to a microscopic theory. In 
section III, using a version of the Dyson method and 
results of the reduced propagator approach we 
show that whenever the environmental evolution is not 
specified and interaction with a driving laser field is 
excluded, the spectral diffusion model is valid, while the 
correlation function {r]tT] T ) is proportional to the real 
part of the environmental correlation function /3r(t — r). 
In section III, where the coupling to the laser is restored, 
we notice that the real obstacle to the rigorous analysis 
beyond the Markovian limit within the phenomenological 
theory is the independent rates of variations assumption, 
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which may be overcome using the microscopic methods 

of [USt, 

as done in section V. Finally, in section VI 
we establish the analytic relations beyond the Marko- 
vian approximation between the single emitter photon 
statistics and the thermal environmental correlation 
function, using a generalization of the generating func- 
tion technique for the photon counting events. Section 
VII summarizes the main steps and concludes the article. 



II GENERALIZED LANGEVIN EQUATION 

In this section we review the theoretical limitations of 
the spectral diffusion approach and compare the Liouvillc 
equations obtained by the phcnomcnological and the mi- 
croscopic approaches. For concreteness, we consider the 
two-level system approximation of a single emitter, which 
is embedded in a thermal environment. For our current 
purpose we exclude the interaction of the emitter with a 
monochromatic laser pump and the electromagnetic vac- 
uum, which will be taken into account later. Thus, the 
unperturbed Hamiltonian of the particle is 



(1) 



where &i denotes the Pauli matrices. Within the spectral 
diffusion approach the two-level system density operator 
p{t\r) t ) is a function of a real process r] t , and obeys a 
stochastic Liouville equation [231-425 1 



— p{t\rj t ) = i 



H S + -jTit& z 



(2) 



Rewriting Eq. ^ in a vector notation denned by p = 
(pec p gg , Pgc p C g) T , where the suffixes e and g refer to 
the excited and the ground states respectively, we have 



where the superoperator 



f 
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p(t\vt), 



(3) 
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p(t\vt),H s 



(4) 



and 



represents the reversible dynamics 

the superoperator T, reflecting the contribution of the 
random part i [p(t\rjt), hvt&z] > is given by 



f = 



/ \ 



i 

\ -i J 



(5) 



Regarding r)t as white, Eq. ((3j) is a standard Langevin 
equation describing evolution of a Brownian particle (26| . 
The well-known generalization of Eq. §3§ for the non- 



Markovian case of colored noise has the form of 11-1 



dt 



p(t\Vt) 



C + r] t T 



p(t\m) 



l(t, T )p( T \VT)dT, 



(6) 

where the memory kernel 7(t,r) is proportional to the 
correlation function (r)tf] T ), as a manifestation of the 
fluctuation-dissipation theorem. Eq. ([6]) shows that 
the dynamics of the open system generally depends on 
its earlier states, because the environment is capable 
of "remembering" the history of the particle's evolu- 
tion. Since r\ t and j(t, r) cx (ffc^r) are related by the 
fluctuation-dissipation theorem, it is inconsistent to 
manipulate the former without altering the latter. Note 
that although in the Markovian limit, {T}tr] T ) oc 5(t — r), 
Eq. (|6|) becomes local in time, Eqs. ([3][5]) cannot be 
attributed to the Markovian limit of Eq. ([6]) because the 
latter generically includes dissipation arising from the 
last term, which is not reflected in Eqs. ([3][5]). Thus, for 
a rigorous treatment of an open system interacting with 
the thermal reservoir beyond the Markovian approxima- 
tion, a priori, we cannot use Eq. ([3]) with rjt colored, but 
must start the analysis from Eq. ([B]). 

The phcnomcnological spectral diffusion approach was 
initiated because the exact Hamiltonian of the entire sys- 
tem is often unknown, making unavailable the projection 
of the time evolution equation of the total system on the 
particle subspace. Yet, many systems may be satisfacto- 
rily described by a simplified model Hamiltonian. A judi- 
cious choice of such a Hamiltonian is capable not only of 
providing the microscopic analogs of Eqs. (|3|6[) . but also 
of shedding light on the microscopic origin of the spectral 
diffusion parameters. For the particle ® reservoir system 
it is reasonable to choose the Lindblad Hamiltonian 17 1 



H= He 



Hi, 



A 



Db x + Lb{ 



(7) 



where L, L' are the particle operators in the Schrodinger 
picture, b\, b x are the annihilation and creation operators 
of bosons in the reservoir mode A, and g\ are the coupling 
coefficients. Projecting the propagator of Eq. on fixed 
initial and final environmental states, using the so-called 
reduced propagator approach 1^- 21 1 , within the second 
order approximation in the coupling strength g\, yields 
the following master equation for the particle: 



dt' 



p(t\z t ) = i p(t\z t ),H s +z* t Lp{t\z t ) 
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-z to L* p(t\z t ) + z t p(t\z t )L^ - z* o p(t\z t )L- 
drpr (t - r) Dc~ i&s ^-^L P (t\ z T )c l " s (*" T ) - 
dT^(t-T)e- i6s ^-^p(T\z T )Lh i6s ^L. (8) 



Here, Zt is a complex function representing a time depen- 
dent average of the environmental states, weighted by g\ 
Til- 2l| (see also Appendix B), and 



/3 T (t)= / dujJ{u) 



coth 



2k b T 



cos(u;£) — i sin(cji) 



(9) 



(in the continuum limit A — > u>) is a transform of the 
environmental spectral function 



J(w) = g(u) 2 D(u), 



(10) 



where D(u>) is the density of the environmental states. 
Integrating over the initial and the final reservoir states 
one can show that I19h21| 



(11) 



(zt)=0, (z t z T )=0, (z t z* T ) = pr(t-T), 



implying that Zt may be interpreted as a random zero 
mean process with a correlation function /3t(*) (Eq. ©). 
Hence, the reduced density matrix p(t\z t ) in Eq. ([8} is 
conditioned by a given environmental evolution path, 
analogous to the phenomcnological density matrix 
p(t\rjt) in Eqs. (|2l6p . conditioned by a realization of the 
spectral diffusion process rj t . In other words, Eq. © is 
the microscopic counterpart of Eq. ([6]) . 

To demonstrate a disagreement between the spectral 
diffusion method Eqs. and the microscopic approach 
Eq. ([S} we shall specify the coupling operator L. Since 
the interaction of the emitter with a thermal environment 
is usually insufficient to generate transitions between the 
states of an unperturbed system, it is manifested as a ran- 
dom noise perturbing only the energy levels. This effect 
accounts for a self adjoint coupling L = U = K, where 
for a two-level system K = a z . With such a substitu- 



tion, well-known as the spin-boson model [171 ] , rewriting 
Eq. (j8J) in vector notation leads to a Liouville equation in 
the form of Eq. ((BJ, where r\t => A Zt = (z t * — z to ), C = £ 
is given by Eq. (@}, 



/Re[A 2 J \ 

-Re[A z J 
-Tm[A 2 J 

\ nm[A gt y 



(12) 



and the memory kernel j(t) is given by 



j(t) = -2Re[/? T (t)] 



/ 1 

1 

e tU)ot 

\ e~ lUot ) 







iujQt 



(13) 



Superposing Eqs. (j3ll5)) and Eqs. (|6I4I12I13[) in the Marko- 
vian limit, we see that the only way to make the phe- 
nomenological and the microscopic methods agree is by 
restricting A Zt and /3y (t — r) to be purely imaginary. In- 
specting Eq. ([§]) we see that the latter constraint may 
be only approximately satisfied for T — > 0, which corre- 
sponds to the non-Markovian regime. 



Ill DYSON EQUATION 

In this section we show that despite the discussed 
limitations of the phenomenological approach described 
above; the spectral diffusion model Eq. ([3]), with rj t white 
or colored, may be justified under standard experimental 
conditions, while (jyt^) may be identified, up to a con- 
stant, with the real part of the microscopic environmental 
correlation function Prit — t). Note that in practice the 
evolution of the environment is usually not determined 
on the microscopic level, so that the quantity measured 
is not p(t\rjt) [or p(t\z t )], but rather its mean value, given 
as 

Ps (t) = (p(t\Vt)) [or p s (t) = (P(WL (14) 

which is the averaged reduced density matrix standardly 
used in the literature. For a general stochastic equation 
in the form of 



j t p{t\Vt) = 6 [p(t\Vt), t] + r,t?p(tfa), 



(15) 



where O is any functional of p(t\r}t), the equation of mo- 
tion for the expectation value ps(t) = (p(t\r) t )) may be 
obtained using an analogue of the Dyson method. We 
expand the propagator of Eq. (|15[) by iteration as 



U(t,t \Vt) = J2u {n) (t,t \Vt), 



(16) 



n=0 



where 



dt, 



to 



t„ rt 2 

dt n -x... / dtix 

t Jtn 



x U (t, t n ) Vtn tu {t n , t n - 1 )r )tn _ 1 t...T hl tUo(t u to), (17) 
and Uo(t,t') obeys Eq. (JTSJ) with T)t = 0, i.e., 



-u (t,t') = d[u (t,t'),t 



(18) 



Each summand in the rhs of Eq. (|T6|) suggests the well- 
known interpretation of the "free" evolution, generated 
by the unperturbed propagator Uo(t,t'), interrupted 
n = 0,1,2,... times by the random "potential" ?7tT. 
Acting with the average of Eq. (|TB|) on the initial state 
vector p(to) yields an expression for ps(t) in terms of a 



4 



sum of integrals including the zeroth, first, second, etc. 
moments of r\t , as shown in Appendix A. 

Further progress is possible if (r]t n r ]t n -i--- r ]ti) may be 
factorized into a product of correlation functions of a 
lower order, as, for example, occurs for a Gaussian noise. 
Due to the special properties of the latter, assuming it 
is of zero mean, all the odd moments vanish, while every 
U( n) (t,t \vt) defined by Eq. dTTJ) with n even gives rise 
to n\j (2^2.!) terms, differing one from another only in 



the way (^t„^t„_ 1 •••ihi) i s factorized [26[. These terms 
may be represented by Feynman diagrams and classified 
according to their physical meaning. For example, 



U {4) (t,t \vt) 



t pt4 pt3 pt2 

dt 4 d t3 d t2 dtAnttVta^ilti}^ 

3 JtQ Jto Jto 



x U (t, U)TU (U, t 3 )TU {t 3 , t 2 )TUo(t2, ti)TU (h,t ) 

(19) 

leads to three different diagrams resulting from 



(20) 



Now, if the magnitude of the correlation function {rjtrj T ) 
is likely to decrease with an increase of t— r, which can be 
true even in the non-Markovian regime, it is clear that the 
first summand (ffeffe) (vt 2 r ]ti) on the rhs of Eq. (|2"0} gives 
rise to a diagram which dominates over the other two di- 
agrams arising from (vuVt^ivtsVt^) and {n] ti Vt 2 ) {Vt s i}ti) > 
because the time interval between the subsequent mo- 
ments is always the smallest. Hence, to a reasonable 
approximation, all the diagrams including the correla- 
tion functions of a pair of non subsequent times may be 
neglected, which leads to a master equation 

p S (t) = 0[ps(t)A+ / M(t,T)p S (T)dT, (21) 



(It 



whose memory kernel is given by 



M(t,T) = (T) t TI T )?U (t,T)r. 



(22) 



This quite general result means that the effect of a 
colored noise, arising from an interaction with an 
environment and rigorously described by Eq. ([6]) would 
be indistinguishable from that predicted by Eq. ((3|), 
provided the latter is constructed such that the resulting 
Dyson equations (f2"Tj) coincide. 



To apply Eqs. (|21I22[) to the spectral diffusion model 
Eqs. (J3]l5j) , we set O [p(t\r]t), t] — Cp(t\r]t), where C is 
given by Eq. (jH). In such a case the solution for Uo(t, t') 
becomes trivial and yields 



/0 



M(t,r)= -{mVr) 







e iW °(*- 
\o 



\ 







-iu (t—r)j 



(23) 



On the other hand, setting L = U = K and assuming the 
reservoir is initially prepared in Boltzmann equilibrium 
at temperature T, eliminating the environmental degrees 
of freedom in Eq. (0) leads to [H-Eil 



TtPs(t) 



H s ,ps(t) 



dr- 



dr, 



( 24 ) 

Substituting H$ = jUJo(t z , K = b z and rewriting 
Eq. (|24p in vector notation we obtain an equation in the 
form of Eq. J5TJ), where 6 [p{t\r) t ),t] = Cp\t\rj t ), with L 
given by Eq. (j4]), and 



M T (t-T) = -4Re[/3 T (t-T)} 



/0 \ 



e iWo( ^ T ) 

yo o o e-^'^-^J 



(25) 

which constitutes the microscopic analog of the phe- 
nomenological memory kernel Eq. (|2"3"]l . We see that in 
case of a stationary process 1 the two coincide, provided 



4Re[0r(t-r)] = (vtVr)- 



(26) 



In such a way, under all mentioned approximations, the 
phenomenological model Eqs. (|3j[5j) , driven by a random 
real spectral diffusion process 77*, white or colored, leads 
to the same marginal master equation for the reduced 
density matrix ps(t) as the microscopic approach. There- 
fore, despite the limitations discussed in the previous sec- 
tion, it may serve as a shorter effective form for the de- 
scription of the two-level system conditional probability 
density matrix dynamics, whenever the environmental 
evolution cannot be fixed on the microscopic level. In 
other words, the spectral diffusion correlation function 
(r]tri T } may be assumed to be some arbitrary function of 
time (i.e., not necessarily a 6- function) , as it may be rig- 
orously identified with the real part of the environmental 
correlation function, i.e., 4Re [/3t(^ — t)]. 



IV INDEPENDENT RATES OF VARIATION 
APPROXIMATION 

Under the assumptions discussed in the previous sec- 
tion we saw that when the environmental evolution is 



[1] The model Eq. J7J is sufficient to yield only a stationary en- 
vironmental noise, whose memory kernel /3x(t — t) Eq. J9]l is 
invariant under time translation. An extension of the method to 
non-stationary noises may be achieved by a "nested doll" model, 
where the particle ® reservoir system is itself considered as a 
subsystem of a larger environment. 
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not determined, the spectral diffusion model yields a 
master equation matching the one obtained by the ex- 
act microscopic analysis. Thus, the inconsistency of 
the phenomenological equation ([3]) with the fluctuation- 
dissipation theorem in case of colored noise is effectively 
eliminated after averaging over all the realizations of r)t, 
and the Dyson equation (|2"Tj) is then valid also beyond the 
Markovian limit. Let us, however, recall that such a con- 
clusion has been obtained assuming Hs oc L = L', i.e., 
excluding single emitter interaction with the driving laser 
field. Restoring the laser pump within the rotating wave 
approximation [lij . the two- level system Hamiltonian is 
no longer diagonal in the eigenbasis of Hs introduced by 
Eq. ([IJ, and is given by 



1 



H s(t) = ~u v z + cos(w L t) [6+ + . 



(27) 



where lul and 



d c „£ 



are the angular and the Rabi 



frequencies of the laser (d Gg = d ge are the matrix ele- 
ments of the off-diagonal electric dipolc moment, and £ 
is the laser amplitude). Switching to a rotating frame by 
the unitary transformation R(t) = exp [i^t&z] (2f|, al- 
lows eliminating the explicit time dependence, and yields 

H s = \^l<J z + ^fio(o+ + o-_), (28) 

where A^ = luq — u>i is the detuning, which entails the 
corresponding Liouville operator 



C = 



( 
























' 2 




\ 


,• Sk. 

1 2 


1 2 



2 

,-n 



* 2 

-if- i&L 



^ 





-i&L 



(29) 



According to the spectral diffusion approach, the master 
equation for the particle is now obtained by Eq. ([3]) 
where C is given by Eq. (|29p. while the interaction with 
the thermal environment, described by r\t and T Eq. ([5]). 
stays unaltered 2~> . It is evident that by simple 
addition of a real random process to the unperturbed 
emitter transition frequency loq, the phenomenological 
model independently imposes the rates of variation 
associated with the interaction with the thermal bath 
and the laser, as if each coupling acted alone. In other 
words, the spectral diffusion approach is based on 
the so-called "independent rates of variation" approx- 
imation, neglecting the effect of possible correlation 
between the laser and the thermal reservoir, induced by 
the coupling to the two-level system. Generally, this 
approximation is legitimate when the time scale of the 
particle evolution induced by the coupling to the laser, 
i.e. fig j is much longer than the correlation time r c 
of the relaxation process (determined by the deca y o f 
(VtWr)), and becomes strictly exact only if r c — > [161 ] . 
This statement is supported by the microscopic analysis, 
since considering the terms e _zffs ^ _r - ) I/p(r|2; T )e lffs ^ _ ' r - ) 



and c~ 4ffs (*~ r )/3(T|z T )I/'!'e li?s ( t ~' r ) in the integrands of 
Eq. (|8"|). it may be readily noted that an incompatibility 
of the particle Hamiltonian Hs with the Lindblad 
operators L, £J has an effect only beyond the Markovian 
limit. Therefore, the independent rates of variation 
assumption, intrinsically inserted into the phenomeno- 
logical spectral diffusion method, constitutes an essential 
consequence of the Markovian approximation, and it is 



natural to expect that for 



H S ,L 



7^ the phenomeno- 



logical Dyson equation (|21[) will no longer coincide with 
the analogous equation (|24[) obtained by the microscopic 
analysis beyond the Markovian limit. 

To sec explicitly how the independent rates of variation 
approximation alters the results of Section III, we first 
revise the phenomenological scheme, outlined in Eqs. (12 14 - 
122]) . Setting 6 \p{t\t] t ), t] = £/?(%), where now t is 
given by Eq. (|29|). we arrive at a master equation in the 
form of Eq. (f2~Tj) , whose memory kernel is 



M(t): 



/0 \ 



B(t) C(t) 

\0 C{t) B*(t)J 



(30) 



where, using the definition of the generalized Rabi fre- 



quency £1 = y/K 
Bit)-- 



2il 2 



+ COS [fit] (fig 



2A|) 



O 



2iA L tt sin [fit]] , 
(31) 



CW = <^o>§sin 2 



On the other hand, substituting Hs, given by Eq. 
into Eq. (|24p (with K = a z ), we arrive at a master equa- 
tion in the form of Eq. (|2ip . where C, once again, is given 
by Eq. (f29|) . whereas the memory kernel takes the form 



M T {t) 



where 



/ \ 



A T (t) A T (t) Br(t) C T (t) 

V A* T (t) A* T {t) C T (t) B* T (t) ) 



fio 



(32) 



A T (t) = 2ilm Writ)] ^ {A L (1 - cos [M]) - iQ sin [M]} . 



B T (t) = - 



2Rcj3 T (*)] 



{^ + (2A| + «o)cos [fii] + 2iA L ^sin [VLt]} 



C T {t) = 4Re[/3 T (t)]^f sin 2 



n 
—t 

2 



(33) 



Remembering that (rjtrjo) = 4Re[/3r(i)], we compare 
Eqs. (15011551) and see: B(t) = B T (t), C(t) = C T (t), 
while the difference enters through Ax(t), which vanishes 
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in the phenomcnological case. Note that Ar(i), unlike 
Bx{t) and Cr(t), is proportional to Im^^)]. In the 
Markovian limit the real part (dissipation) of the envi- 
ronmental correlation function dominates over the imag- 
inary part (fluctuation). This means that Ax(t), induc- 
ing disagreement between the microscopic theory and the 
spectral diffusion model, becomes important in the non- 
Markovian regime. The Markovian limit of Eqs. (|30|32j) 
may be recovered by setting t = 0, which yields 



M(0) = M T (0) = -4Re [(3 T (0)} 



/ 0\ 



10 

\ 1 J 



Hence, in the Markovian approximation the phenomeno- 
logical method gives a result which coincides with the 
microscopic approach, as expected. In summary so far, 
we conclude that due to the independent rates of vari- 
ation assumption (and not because of the apparent in- 
compatibility with the fluctuation-dissipation theorem) 
the spectral diffusion model Eq. ((3j) does not consistently 
describe the SMS experiments beyond the Markovian ap- 
proximation. The microscopic approach, on the other 
hand, permits avoiding the independent rates of variation 
assumption, and also clarifies the role of the environmen- 
tal spectral function J(w) Eq. (fit))) and the temperature 
of the thermal reservoir T. 



V MASTER EQUATION FOR A TWO-LEVEL 
SYSTEM INTERACTING WITH TWO 
RESERVOIRS 

In this section we set up a reduced master equation 
describing the SMS experiments beyond the Markovian 
limit using the microscopic methods of [Bill. Recall 
that we considered a two-level system emitter contin- 
uously driven by a classical monochromatic laser field. 
Were complete isolation from the environment realistic, 
the particle would undergo simple Rabi oscillations 
[B|. In practice, the unavoidable interaction with 
the electromagnetic vacuum induces incoherent decay 
transitions accompanied by events of the spontaneous 
photon emission, while the coupling to the thermal 
environment is associated with the origin of the spectral 
diffusion noise. 



The total system, in the rotating waves approximation, 
is described by the Hamiltonian (see Appendix B) 

F=^a z + ^(a + +a_)+E M K-^)4^+Ex^SiSxH 

+ E P Pn (o+ «m + °- 4 ) + E a 9\°z (b{ + b x ) , 

, (34) 

where the first two terms describe the Hamiltonian of 
the two-level system driven by the laser; the 3-4th terms 



are the free Hamiltonians of the vacuum and the thermal 
bath, whose interaction with the particle is described by 
5th and the 6th terms respectively. The theory of the re- 
duced propagator and its application to the derivation of 
the marginal master equation was comprehensively de- 
veloped for a particle interacting with a single bosonic 
field in 

(Bill- 

In case of a particle simultaneously cou- 



pled to several baths, the derivation of the master equa- 
tion might be complicated by the entanglement arising 
between the reservoirs from the induced indirect interac- 
tion. However, this interaction, as shown in Appendix 
B, is manifested in terms going beyond the second order 
in coupling strength, and may be neglected if the inter- 
action is weak. Assuming in our case that the two-level 
system is weakly coupled to both reservoirs, it is justified 
to approximate the reduced evolution of the particle im- 
posing the contributions arising from the interaction with 
each one of the reservoirs independently, which yields 

d ;Ps{t) = i \ps(t),H s - 



dt' 

— / dra(t — t) 



<r+,e 



-iH s {t-r). 



JH s (t-r) 



dra*(t-T) 
dT/3 T (t - t) 



r ^s( t -r )/5s(T) ^ +e ^ s (t-r) ; ^ 



-iH S (t- 



dr^it.-T) 



iH s (t-T. 



As(t)<7 2 



iH s (t- 



(35) 

where Hs is given by Eq. Taking the Markovian 

limit of the electromagnetic vacuum correlation function 
a{t — t) (x 5(t — t) and rewriting Eq. (|35[) in vector 
notation, we find 



d , . 



Ps{t)+ / M T (t-T)p s (T)dT, (36) 



In 



where 



C-- 



(-1 


1 2 



,-fio 



iA, 







oho 



(37) 



r 



o\ 
o 



0/ 



(38) 



/ 
7 

\0 

7 is the spontaneous emission rate B> 3 , and the ther- 
mal memory kernel A4t(£), already calculated, is given 
by Eqs. (|32I33|) . Eq. (|36| is the desired reduced master 
equation, unrestricted by the assumption of independent 
rates of variations, for a driven two-level system interact- 
ing with two bosonic reservoirs: one within the Marko- 
vian approximation (the electromagnetic vacuum), and 
another beyond it (the thermal bath). 
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VI ENVIRONMENTAL CORRELATION 
FUNCTION VS PHOTON STATISTICS 



In what follows we employ Eq. (|36|) (together with 
Eqs. (|32l37l38p ) for the investigation of the thermal en- 
vironmental correlation function using data of the single 
molecule photon statistics. Several methods for count- 
ing spontaneous photon emission events were proposed 
in the past. Prominent examples are the quantum jumps 
approach 27 1. or the generating function method [23J, 
which combines analyticity, calculational simplicity (for 
the low-dimensional systems) and intuition. Even though 
originally developed for the Markovian Optical Bloch 
Equations [lil ] , the generating function approach may be 
quite easily generalized to the time-retarded equation of 
the type of Eq. ([3T>)) pjj . Setting to = and transforming 
Eq. (f3"6")l to Laplace space t — > £ we have 

Ps(0 = [C-£- M T (0 - f ] _1 Ps(0) = G(OPs(0). 

(39) 

The propagator G(C) may be formally expressed in terms 
of the iterative expansion 

G(0 = [i + &(C)f + Go(C)m(C)r Go(C) = 



where 



and 



l-&(C)f] Go(C) = J2G in) (C), (40) 

n=0 

G {n) (C) = (Go(Or) n Go(C), (41) 

Go(0 = k-C-MriOV 1 (42) 



is the propagator of Eq. (|39[) without f . The operator 
T, given by Eq. (|3"5)) . couples the population of the 
excited state p cc directly to the population p gg of 
the ground state, and hence, describes an incoherent 
transition, which (within the first order perturbation 
approximation with respect to the linear coupling to the 
radiation field) may be associated with the spontaneous 
emission of a single photon [29|, [3(| ■ With such an 
interpretation, each term "'((,') of Eq. (140!) corresponds 
to the two-level system evolution conditioned by n 
spontaneous photon emission events. 

The iterative expansion of Eq. (|40[) may be compactly 
represented using the generating function (23j 



G(t, 



n=Q 



(<), 



(43) 



where 

p< n Ht) = G {n) (t,0)m, (44) 
and the parameter s serves as an odometer of the spon- 
taneous emission events. Substituting the definition 
Eq. (HU into Eq. @g) yields 



G(C, s)=[(-C- M T (C) - st\ G(0, s), (45) 

which differs from Eq. p9p only by the extra factor s 
multiplying T. Since 



s=0) 



(46) 



the statistics of the photon emission events may be fully 
obtained by the derivatives of G(t,s) with respect to s. 
The most common and simply measured quantities are 
the probability of n spontaneous photon emissions 



^(0=^(0+^(0= 



d n - d n - 

^ G ee( S ,C) + ^G gg ( S ,C) 



(47) 



s=0 



the mean photon number 

oo 

and the second moment 



dG cc (s,0 , dG eg (s,() 



ds 



ds 



8=1 

(48) 



(^(C)> 



a 2 G ee ( S ,C) , d 2 G gg ( S ,0 



ds' 2 



ds 2 



dG ec ( s ,C) , dG gg (s,C) 



ds 



(49) 



Using Eqs. (|37I38I32I33I45I47I48I49|> it is now straight- 
forward to determine the connection between the thermal 
noise correlation function Pxit) and single molecule pho- 
ton statistics. For this purpose we first need to find the 
Laplace transform of the thermal memory kernel A4x(t) 
Eqs. (|32l33p . It follows from Eqs. ((33J) that the matrix 
elements of M.r{t) may be decomposed as 
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A T (t) = i& 21m\flr(t)]$£ + c' Ml lm[p T (t)] 



_£2 Ax 



IM*)] (& 



1 + 2 



BTW = -i[2Re[^(t)]+e-^Re[^(t)](l+2^(^-^))+e^Re[/3 r (f)] r , ^ ^ , ^ 

C T (t) = § [2Re[/3 T (t)] - e- <tn Re[/? T (i)] - c* to Re[/3 T (t)]] . 
Since for a function of time /(f), we have L t ^Q [e ±lf2 '/W] = /(C + in), the Laplace transform of Eqs. (|50]) yields 



(50) 



MO 



'n 2 



2^Im[/3 T (C)] + Im[/3 T (C + - - Im[/3 T (C - + ^) 



5t(C) = 



o- 

"IP 



2Rc 



^ T (C)]+Re[M+^)] (l+2^(^-s|))+Re[/MC- in)](l+2^(^ 



(51) 



Cr(C) = 



o 2 



2Re[^ T (C)]-Re[/3 T (C+^)]-Re[/3 T (C- i O)] 



The explicit calculations following the prescription Eqs. (|37I38I32I45I47I48I49I51|) may be done with the help of 
computational programs such as Mathematica. To illustrate the method we confine the following example to frequently 
used experimental conditions. Assuming that at t — the emitter is prepared in the pure ground state and the laser 
frequency is close to the resonance; expanding the results in a Taylor series in A/, up to the first order wc find 



G cc (C,s)+G gg (C,s) 



3a + 7 + 2C) (C + 7)+202 



+A, 



C (85 + 7 + 20 (C + 7) + (7(1 - «) + 2C) n 2 
J(s-l)7((7+C)(8o+7+2C)+2ng) (26(7+20+8o + (6-6_)-(7+2C+2in )&-+8a_(6-6+)-(7+2C-2ino)5 + ) 



(C( 7 + 0(85 + 7 + 20 + (7(1 - s) + 2C)figr (7 + 2C + 4(5_ + 5+)) 



where we used the shorthand notation: a = 5(0 = 
Re[/3(0], b = 6(C) = Im[/3(C)], a± = o±(C) = MP(C ± 
ifi)] and b± = 6±(0 = Im[/3(C ± in)]. Eq. {52]) indicates 
that on resonance A^ = the photon statistics depends 



on a 



Re 



/3(C) 



alone, which makes the reconstruction 
of the real part of the thermal correlation function partic- 
ularly simple. Substituting Eq. (J52J into Eqs. (|47l48l49j) 
in case of on resonance excitation yields 



Pn(C) 



(n(0) 



n! 7 "ng" (( 7 + 0(7 + 2C + 85) + 2^) 

(( 7 + 20( 7 C + C 2 + n 2 ) + 85C( 7 + Of +1 ' 
7 n 2 



C 2 ((7 + 0(7 + 2C + 85) + 2ng)^ 



( - 2 7^(7 + C)(7C + 2(C 2 + ni) + 8a^ 

C 3 ((7 + 0(7 + 2C + 85) + 2n2) 2 

On specifying J(ui) and the temperature T, which ac- 
cording to Eq. © are needed for the calculation of the 
real part of the thermal environmental correlation func- 



tion Re [Pr(t)] = Jq 00 doj J(uj) coth 
straightforward to use Eqs. ([5j 



cos [ojt] , it is 



2k b T 

for predicting the cor- 
responding photon statistics. The converse procedure of 



reconstruction of the thermal correlation function is ob- 
tained by inverting Eqs. (|53[) . For example, expressing 
5 = Re[/3(0] in terms of Po(C) an d ("-(0): which are un- 
derstood to be measured experimentally, with the help 
of Eqs. (f53|) respectively yields 



Re[8(C)] = 



(7+0(7+20- 



-2ng- 



•(7+2C)(C(7+C)+ng)P (C) 



8(7 + C)(Cfo(C)-l) 



(54) 



Re[/3(C)] 



7 ng-(n(0)C 2 ((7 + 0(7 + 20 + 2ng) 
8(n(0)C 2 (7 + 



(55) 

Note that Eqs. ()54|55|) establish a well-defined relation 
between Po(C) an d (ft(C))- Expressing 5 = Re[/3(0] 
in terms of P n (C) with n > and (5 2 (0) is slightly 
more complicated and may entail multiple possibilities, 
which must be examined from the physical point of 
view. These expressions, which we shall skip since they 
are massive, may be easily obtained with the help of 
Mathematica if needed. 

Finally, we would like to illustrate the results of 
Eqs. (1551) . To do so, using Eq. ©, we must choose 
an explicit form of the environmental spectral function 
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0.02 0.04 0.06 0.08 0.10 0.12 0.14 



FIG. 1: (Color online) The plot of Pi(C) (Eq. (g3|) for 
7 = 1, fio = 2, corresponding to environmental correlation 
function /3(f) = u 2 e~ Rt for v = 1 and i? = 0.1,1,10,100 
represented respectively by the red, blue, green, and orange 
curves. 




of 



sin[Bf])e" 



FIG. 2: (Color online) (a) The plot 
Ae~ m ; Ae~ m2 ; (1 - A cos[Bt])e~ Rt ; (A - 
for A = 0.5 and B = 10, represented respectively by the 
red, blue, green, and orange curves; (b) The plot of -Po(C) 
(Eq. (|53|) h for 7 = 1, fio = 2, resulting from a different choice 
of the environmental correlation function in Fig. 2(b). The 
colors of the curves are matched . 



J(oj). The latter, if microscopically unknown, may be 
modeled phenomenologically, for example, as a power-law 
J{lo) = ?7 s a/'cj^~ s c~"/" c , where r] s is the viscosity coeffi- 
cient, uj c is a cutoff frequency and s, for the case of the 
spin-boson model, is the dimension of space [l?}- How- 
ever, the integral of Eq. © induced by such a choice of 
spectral function, is not analytically solvable in paramet- 
ric form. Hence, we restrict the illustration of Eqs. (|53|) to 
some arbitrary examples of the environmental correlation 
function. In Fig. [1] we plot the probability of one photon 
emission, setting j3(t) = v 2 e~ Rt , which is a standard as- 
sumption of the spectral diffusion approach [23, . The 
graph, comparing A(C) f° r R = 0.1, 1, 10, 100 (red, blue, 
green, orange) shows that an increase in the order of R 
requires nearly exponential improvement of the measure- 
ment accuracy. In Fig. [2] we examine the dependence 
of the probability of zero photon emission events -Po(C) 
(Fig. Wip)) on several shapes of /3(f) (Fig. (2Ua)), chosen 
such that the autocorrelation time r c and the maximal 
value of (5(t) are close. The fact that for different choices 
of /3(t) we arrive at effectively the same result for -Po(C) 
may be seen as a justification of the phenomenological 
assumption /3(t) = v 2 e~ Rt , used to plot Fig.Q] 



VII SUMMARY 

By superposing an approximate phenomenological 
Dyson equation with a microscopic reduced master equa- 
tion, we provided a correspondence between the micro- 
scopic and the spectral diffusion approaches. It was 
shown that excluding the interaction of the two-level sys- 
tem emitter with the laser, the dynamics governed by 
the spectral diffusion model is indistinguishable from the 
dynamics obtained by the microscopic theory, whenever 
the evolution of the environment is not determined on 
the microscopic level. This allowed to identify, up to a 
constant, the spectral diffusion correlation function with 



the real part of the thermal environmental correlation 
function. Furthermore, it was demonstrated that the 
real problem of using the spectral diffusion method for 
the description of standard SMS experiments beyond the 
Markovian limit, is the independent rates of variation as- 
sumption, which can be overcome using the microscopic 
reduced propagator theory Finally, combining 

the methods of the latter with the generalized generat- 
ing functions approach, we have established analytic ex- 
pressions for the thermal correlation function allowing to 
extract information on the environmental noise beyond 
the Markovian approximation from the measured single 
emitter photon statistics. 



Appendix A: Derivation of Dyson equation 

In this appendix we provide the details of derivation 
of Eqs. (1211221) . Substituting p{t\r) t ) = U(t,t \r]t)p(t ), 
where U(t,to\r)t) is the propagator, into Eq. ()15[) we have 



-U{t,t \ Vt ) = O U(t,t \Vt),t +VtTU(t,t \ Vt ), (AT) 
at l i 

and rewrite the latter equation in the integral form fl6| 



U(t,t \r) t ) = U (t,t )+ / U (t,T)ri r tU(t,t \m)dT, 



(A.2) 

where C7q (t, t ) is the solution of Eq. (|18p . Further, we 
expand Eq. (|A.2[) by iteration 



U(t,t \r) t ) = U (t,to)+ / dtiU (t,t 1 )T hl TU (t 1 ,t Q )+ 
+ / dh dt 1 U {t,t 2 )7 1 t 2 fU {t 2 ,h)r hl fUo(t 1 ,t Q ) + 
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+ / dt 3 f 3 dt 2 fdt 1 U (t,t 3 )r] t3 tU (t 3 ,t 2 )ri t2 'Tx 

Jto J to J to 

xu (t2,ti)m 1 ruo(t 1 ,t Q )+ 

J to J to •* to J to 

x Vt3 rU Q (t 3 ,t 2 ) Vt2 fu (t 2 ,t 1 )i hl fu Q (t 1 ,t Q ) + ... . (A.3) 

Assuming rj t is a zero mean Gaussian noise, we take an 
average of Eq. (|A.3[) . Discarding the long time correla- 
tions contribution, as discussed in Section III, this gives 

U(t,to) = U (t,t )+ 

t pt 2 

dh / dh(r]t 2 Vt 1 )Uo(t,t 2 )rU {t 2 ,t 1 )tUo(h,to)+ 

to J to 

ft rt 4 rt 3 pt 2 

+ \diA dt 3 dt 2 dt 1 (rjt 4 'i]t 3 }(rit 2 'nt 1 )Uo(t 1 t4)x 

J to J to J to J to 

( A - 4 ) 

where U(t,t ) = (U(t,to\r]t))- Evidently, the above iter- 
ative expansion is equivalent to 

U(t,to) = U (t,t )+ 

ft pt 2 

+ dh «M%%>#o(«,t2)T0b(fc,ti)T#(ti,to), 

Jt„ J t 

^-r^ (A : 5) 

which in turn, by analogy with Eqs. (|A.1IA.2[) . is equiv- 
alent to the differential equation 

-U(t, t ) = U(t, t ), t + / dT(r) t r) T )rU (t, t T )TU{t T , t„). 

(A.6) 

Finally, acting with the last equation on the initial state 
p(t ) yields Eqs. (|21I22|) . 



where b\ , b\ and a\ , d\ are the boson ladder operators of 
the thermal environment and the electromagnetic field 
respectively, the Pauli matrices <Xj represent the particle 
operators, and f2o is t he Rabi frequency of the laser pump 
oscillating at frequency ojl- Transforming to the rotating 
frame by 



R(t) = cxp 



(B.2) 



we suppress the time dependence in Eq. (|B.1[) . which 
yields Eq. (fM|) of the article. 



Given Eq. (f34|) . extending the methods of lH-21| for 
simultaneous interaction of the particle with two bosonic 
reservoirs, we consider the total propagator of the system 
in the partial representation picture with respect to Hg + 



(B.3) 

and define the reduced propagator 

G(yt^yt o zt \tt ) = (ytzt\Ui{t,t )\yt o z to ) : (B.4) 

where \z t ) = Y[\ \ z t,\) describes the thermal field in the 
Bargmann coherent states representation (and similarly, 
\yi) represents the state of the electromagnetic field). 
Applying G(y% z^y to z to \Uo) to the initial state of the 
two-level system, propagates the latter such that the 
final state is simultaneously conditioned by specific 
evolution trajectories of both reservoirs. In what follows 
we show that under certain circumstances, the particle 
conditional density matrix pi(ifiz%y to zt \tto) depends 
on the environmental degrees of freedom through the 
stochastic processes yt , Zt and their autocorrelation 
functions a(t),f3(t), defined below. This constitutes a 
simplification of the general case, where higher order 
cross-correlation functions of yt and zt can be involved 
as well. 



Appendix B: Reduced propagators approach with 
two reservoirs 



In this appendix we consider a closed system, com- 
posed of a two-level particle interacting with a thermal 
reservoir, a monochromatic laser field and the electro- 
magnetic vacuum. The total Hamiltonian within the ro- 
tating wave approximation (RWA) with respect to the 
coupling to the electromagnetic field is [l6, 18 - 21 [ 



H(t) = fa z + (a+^ e -i"Lt + a_^c luJLt ) + 

+ Ea u\b{b x + E p w M &t a M + 
+ Ea 9^4b{ + 6a) + E M P^+a^ + v-al), 



(B.l) 



The time evolution equation for G(y£ z£y to z to \tto) is 
obtained from the projection of the Schrodingcr equation 
for Ui(t,t ), given by Eq. (|B~3)) : 



m^(yt^yto^to\ tt o) = (ytzt\§- t Ui(t,t )\y to zt ) 



-i Hs + ^Exdxe^z 



t,A" 



-io" z EA3AC" ? " At (27tz t |6AC//(t,to)|yt ^o)- 

-^+E tl P^ ?;( ^"" L) *(y* f *l«M^( i ^o)|yto^o), 



G{y* t zty ta z to \tt ) 



(B.5) 
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where the last two terms constitute an obstacle for get- 
ting a closed equation for G(y£ z^y to z to \tt ). We use the 
same scheme as above to represent these terms as func- 
tions of G(ifi: z£yt o z to \tt ). Starting with the thermal 
reservoir we have 

(ytZt\b\Ui(t,t )\y to z to } = 
= {ytz^U^t^UY^t^hUiitMMoZt,) = (B.6) 
= (ytZt\Ui(t,to)b\(t,t )\yt z to ), 

where 

h(t,t ) = ^(Mq^aCMMo). (B.7) 
Rewriting Eq. (|B.7|) in the integral form yields 

k(t, t ) =b x - igx I <Mt, t )c^ r dr, (B.8) 

where <j z (t, to) = U~j~ (t,to)a z Ui(t,to)- Repeating the 
same procedure with respect to the radiation field and 
inserting the results back into Eq. (|B.5[) gives 

§-tG(ytz^yt z to \tto) = -iH s + a z z* t - a z z to + 

+^~Vt - p+yt ] G(yl^yt zto\tto)~ 

-io z f to Pit ~ T ){ytZ t \Ui{t 7 to)a z (T, t )\y to z to )dT- 

-iv+ft a{t - T)(y t z t \Ui{t,to)a-{T, to)\y to z to )dT, 

(B.9) 



where 



z T = -i^2g\c luJxt z T ^ x , y T = -i^p^-" ~ UJ "y Ttll , 



and 



0(f) = £ | flA | 2 e- 4 



, a(t)^J2\P^ c ~ 1 ^ 



(B.10) 



(B.H) 



Further, in order to close Eq. (|B.9|) . we are interested 
in rearranging the last two terms in Eq. (|B.9() as 



{y t z t \Ui{t 1 to)a l {T, t )\yt zt ) = (ytZt\o-i(r,t)Ui(t,to)\yt z to ) 
= 0{y;z*y to z to tT)G(y^z^y to z to \tto), 

(B.12) 

where 0(y^ z^y to z to tT) is an operator acting in the par- 
ticle subspace. This can be done by formally integrat- 
ing and iterating the Heisenberg equation for (r, t) in 
power series of gx and p^. Taking into account that the 
autocorrelation functions a(t), j3{t) are already of the sec- 
ond order in gx and p M respectively, keeping only the 
zeroth order solution for <7i(r, i): 



accounts for a satisfactory approximation in the case gx 
and p M are small. Substituting the truncated expansion 
Eq. (|B.13[) back into Eq. (|B.9|) we neglect all the terms 
proportional to g^P™ with m+n > 2. Since the zeroth or- 
der approximation Eq. (|B.13[) is restricted to the particle 
subspace, we can pull it out of the brackets in Eq. (|B.12[) . 
which allows to close the equation for G(y% z£y to z to \tto) 
and yields 

§i G( S*t^yt Q z to \tto) = 



-iHs + o z z* t - o z z ta + a-y* t - a + y t(i 



-ia z J* o 0(t - T ) c -*Hs(t-T) &zel H s (t~r)_ (B.14) 



-ia+ J* a(t - T ) e -iH s (t-T)fr_ e iX s (t-T) 



*G{ylz$y to z to \tto)dT. 



The lesson of Eq. (|B.14[) is that within the second order 
approximation in the interaction magnitudes, the effec- 
tive coupling between the two reservoirs is eliminated. 
The reduced evolution of the particle is given by direct 
independent summation of the contribution arising from 
the interaction with each one of the fields alone. 



From now on the results of may be adopted 

directly. Using Eq. (|B. 14|) the master equation for the 
conditional density matrix follows from 



^tPi(ylzYyt zt \tto) = §i G(y^z^y to Zt \Uo)x 



x Pi($ % ytoZt \toto)G($%yt zt \tt ) 



(B.15) 



Afterwards, the evolution equation for the marginal 
density matrix pg (i) is obtained by integrating out the 
irrelevant environmental degrees of freedom, taking into 
account the initial Boltzmann equilibrium distribution of 
the thermal bath and the zero temperature ^-correlated 
distribution of the electromagnetic vacuum. This proce- 
dure, requiring an application of the generalized Novikov 
theorem, was worked out in (l9l - [2l| . Using the results for 
each of the two reservoirs respectively we obtain Eq. (|35|) . 



a z (T,t) = e - iH s{t-r) a . e iH s {t-r) + ^ 



(B.13) 
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